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SUMMARY 


Based on linearized equations of motion utilizing only the three 
moment equations and assuming only flat-spin conditions, it appears that 
contemporary designs fwith the moment of inertia about the wing axis Iy 


considerably greater than the moment of inertia about the fuselage axis 

(rolling -moment coefficient due to 

(rolling -moment coefficient due to 


1^ having positive values of C^ 
rolling) or positive values of C^ 


sideslip) will probably not have a stable spin in the flat-spin region 
near an angle of attack of 90°. If the damping in pitch in flat-spin 
attitudes is zero, stable flat-spin conditions may not be possible on 
an airplane having the mass primarily distributed along the wings. The 
effect of moving ailerons with the spin or the effect of applying a 
positive pitching moment producing recovery for contemporary fighter 
designs will be greatest for large negative values of C n (yawing- 


moment coefficient due to sideslip). In addition, for a certain critical 
value of positive C R p, the rolling moment applied by moving ailerons 


with the spin or the application of a positive pitching moment will have 
no effect on reducing the spin rate. 


INTRODUCTION 


Flat spins of airplanes have become more prevalent than In the past, 
along with the trends toward lengthened fuselage forebodies, Increased 
relative distribution of mass in the fuselage, and low-span wings. 
Research results. In general, have indicated the important effects of 
such factors as mass distribution and fuselage-nose cross-sectional 
shape on the overall stability of a potential spinning motion in deter- 
mining the nature of spins achieved for current configurations. 
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As regards the effect of fuselage-ncse cross-sectional shape, results 
of investigations of dynamic models in the Langley 20-foot free-spinning 
tunnel and of static tests in the Langley 300-MFH 7- by 10-foot tunnel 
(refs. 1 and 2) have indicated that certain fuselage cross-sectional 
shapes can provide propelling moments at flat-spin attitudes which might 
make an airplane have an uncontrollable flat spin with a high rotational 
rate. Spin-model tests have also shown, however, that, for certain models 
of contemporary fighters which had propelling noses, the flat spin was 
not a stable condition but that the motion of the model appeared to be 
oscillatorily divergent in roll. This oscillatory roll divergence, when 
it occurs, generally causes the spin rate of the model to decrease and 
the model either ceases to spin or assumes a steeper spin attitude with 
a slower rotational rate. The mass distribution of models for which 
this motion has been observed is such that most of the mass is extended 
along the fuselage (the moment of inertia about the wing axis 1^ is 

high ) and the wings are relatively light (the moment of inertia about 
the fuselage axis 1^ is low) . From the results of the aforementioned 

spin-model tests, it is considered desirable that a given airplane be 
unstable in the flat-spin region and that the aerodynamic or mass char- 
acteristics should be arranged so that there is no possibility of a 
stable flat spin. 

Accordingly, an analytical investigation was undertaken to determine 
how a flat-spin condition, on a specific design, would be influenced by 
stability characteristics. This problem is analyzed herein by utilizing 
only the linearized pitch, roll, and yaw equations of motion. This 
approach seems reasonable inasmuch as model force-test data show little 
change in the aerodynamic forces for small sideslip angles and for varia- 
tions in angle of attack in the flat-spin region. Also, as stated in 
reference 3* the recovery motion of the airplane appears to be affected 
primarily by the action of the moments rather than of the forces. The 
effects of applying various moments in a flat spin have been examined. 


SYMBOLS 


The body system of axes is used. This system of axes, related 
angles, and positive directions of corresponding forces and moments are 
illustrated in figure 1. The symbols are defined as follows: 


M x 

^V 2 Sb 

y 


rolling -moment coefficient. 



I£ *1* 


3 


My 

C m pitching-moment coefficient , 

ipV 2 Sb 

2 


c 


n 


yawing-moment coefficient , 


1 2 
=pV Sb 
2 


M X 

My 

M Z 

A 

t 

S 

b 

P 

V 

a 

7\ 


coefficients of characteristic equation (eq. 
rolling moment acting about X body axis, ft-lb 
pitching moment acting about Y body axis, ft-lb 
yawing moment acting about Z body axis, ft-lb 
difference between two values 
time, sec 
wing area, sq ft 
wing span, ft 
air density, slugs/cu ft 
resultant linear velocity, ft/sec 
resultant angular velocity, radians/sec 
roots of characteristic equation 


(A19)) 


ft-lb /radian 


% 


ft -lb /radian 




ft-lb /radian 



(3My i 

— — f t -lb -sec /radian 



4 



= ft-lb-sec /radian 


X,Y,Z longitudinal, lateral, and vertical body axes, respectively 

p,q,r components of angular velocity about X, Y, and Z body 

axes, respectively, radians/sec 

r Q value of r before disturbance, radians /sec 

^X'^Y'3'Z moments of inertia about X, Y, and Z body axes, 

respectively, slug-ft^ 

0 angular displacement of X body axis from horizontal plane 

measured in vertical plane, positive when airplane nose 
is above horizontal plane, deg or radians 

0 total angular movement of Y body axis from horizontal 

plane measured in YZ body plane, positive when clockwise 
as viewed from rear of airplane (if X body axis is 
vertical, 0 is measured from a reference position in 
horizontal plane), deg or radians 


a angle of attack, angle between relative wind V projected 

into XZ-plane of symmetry end X body axis, positive when 
relative wind comes from below XY body plane, deg or 
radians 


P angle of sideslip, angle between relative wind V and 

projection of relative wine, on XZ-plane, positive when 
relative wind comes from right of plane of symmetry, 
deg or radians 


horizontal component of tota2. angular deflection of X body 
axis from reference position in horizontal plane, positive 
when clockwise as viewed from vertically above airplane, 
radians 
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F(a) 

f(d) 




2V 


ac 


n 


'n~ - 


a — 

2V 


characteristic (or stability) equation 

numerator of equation for motion after various forcing 
functions are applied 


a ^ 


ac, 


m 


-'nu 


a 


2V 


acj 

ap" 


ac, 


m 


“e ae 

ac 


<V = 


n 


per radian 


per radian 


per radian 


constant value in A0 forcing-function equation (eq. ( 5 )) 
constant value in Ar forcing-function equation (eq. ( 4 )) 


p ap 


K Ar 

D = A 
dt 


Subscripts and superscripts: 
o initial 

t time 


A dot over a symbol represents the derivative with respect to time 
for example y $ - 
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METHODS AND APPROXIMATIONS 


Three degrees of freedom were utilized in the present investiga- 
tion. The moment equations of motion were used to correspond with these 
degrees of freedom and are presented in appendix A as equations (Al) 
to (A3) . As indicated in these equations, the body axes and principal 
axes are assumed to coincide, and the engine gyroscopic moments are 
neglected. Also, it was assumed that the airplane was spinning about 
its center of gravity in a completely steady spin; that is, p, q, 
and r were zero. For these steady spins, the sideslip angle was 
assumed to be zero and in order to achieve equilibrium at zero sideslip 
a small amount of rolling moment was applied to oppose the rolling 

moment generated by 


p and a small amouit of yawing moment ZNMr? to 


oppose Mg r . 


The slope of the pitching -moment curve M v is the same as the 

x 0 

pitching moment plotted against angle of attack. The plot of pitching 
moment against 9 would have the form similar to that indicated in 
figure 2. In addition, the pitching-moment coefficients were nondimen- 
sionalized with respect to the wing span b instead of the chord. 


The aerodynamic derivatives used were assumed to be constant for 
the flat-spin attitudes assumed, and such derivatives as M Y and M 

A r 


s 


were considered to be small enough to be neglected. 


Because of the fundamental assumptions made for the flat-spin con- 
dition assumed, which would be near an angle of attack of 90 °, the fol- 
lowing conditions existed: 

(a) The rate of descent V remained constant and was the resultant 
velocity (determined by equating the drag at the spin angle of attack 

to the weight of the airplane). 

(b) The rate of yawing r was equal to the resultant rotational 
rate about the spin axis ft or ty. 

(c) The sideslip angle 0 was the same as the wing-tilt angle about 
the body axis $. 

(d) The angle of attack a was equivalent to 90° + 0, with 0 
being the inclination of the longitudinal axis of the airplane above 
the horizon. (Thus, for very flat spins, 0 becomes a small negative 
angle . ) 
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With these conditions, and from expressions of Euler T s attitude 
angles in terms of the angular velocities about the body axes (appen- 
dix A), expressions for p, q, and r in terms of AP, A0, and Ar 
were obtained, with second-order terms being considered small enough 
to be neglected, and are presented as equations (A4) to (A 6 ) . From 
these equations, the equations of motion (eqs. (Al) to (A3)) about the 
three body axes are rewritten as equations (Al 6 ) to (Al 8 ). This pro- 
cedure is used in order that only three variables describing the motion 
remain in the equations of motion, whereas originally there were five. 

When equations (Al 6 ) to (Al 8 ) are placed in determinant form, the 
characteristic equation (eq. (A19)) can be determined and is of the form 

CjD 5 + CgD 1 * + CjD 5 + C 4 D 2 + C^D + C 6 = 0. 

The expression for the various coefficients of the quintic are contained 
in appendix A. The conditions necessary for the roots of equation (A19) 
to be stable (according to ref. 4) are that the coefficients of the 
equation all be positive and that certain functions of the coefficients 
known as Routh's discriminant also be positive. Routh's discriminant 
for the quintic equation is as follows: 

C 2 C 3 - CJC 4 ( 1 ) 




- C 2 ^C 2 C 3 - 

C 2 C 3 - C A) - C 2( C 2 C 5 - ° 1 C 6 )] - C 6( C 2 C 3 - 


( 2 ) 

C 1 C 4) 2 

(3) 


An alternate form of Routh’s discriminant is given in reference 5* 
In order to investigate some of the factors that might affect the sta- 
bility of the flat spin, the stability derivatives were varied in the 
stability equations for an airplane representative of current fighter 
designs as regards mass and dimensions (table I). Various combinations 
of derivatives tried are presented in table II, and the coefficients of 
the quintic equations for these combinations are shown in table III in 


terms of C 


n p 


and C 




Stability diagrams in the form of C n ^ plotted 


P 


against Cj are given in figures 3 to 9* By following the usual pro- 
cedures in stability-analysis work, the constant term in the characteris- 
tic equation (the coefficient Cgj equated to zero is plotted as the 


divergence boundary (where one or more of the real roots becomes posi- 
tive); also, the limiting condition of the three forms of Routh’s 
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discriminant , shown in equation ( 3 ), equated to zero is plotted as the 
oscillatory divergent boundary (real root negative, but a pair of the 
complex roots positive for all cases investigated) . When the constant 
term in the quintic ^the coefficient becomes negative, the motion 

is divergent; and when the limiting form of Routh's discriminant becomes 
negative, the motion has oscillatory instability. 

The technique utilized in the solution of motions after application 
of forcing functions is explained in appendix A. 


RESULTS AND DISCUSSION 


Stability Diagrams 


Cases A and B were computed to show the effect of different rota- 
tional rates on the stability boundaries. The rates chosen for cases A 
and B are similar to those that could be obtained in flat spins, where 
the rate chosen for case A ^r Q = 1.9 radians /sec^ would be somewhat slower 

than the average spin and the rate chosen fcr case B ^r Q = 3*14 radians /sec^ 

would be about normal. For these two cases, plotted as figures 3 and 4 
respectively, the damping in roll C 7 was chosen as a reasonable value 

P 

for a swept-wing airplane based on the oscillation tests presented in ref- 
erence 6 . From previous experience and from oscillation tests presented 

in reference 6 , the values chosen for C m , C^ , and C n are considered 

q 8 r 

to be within the range obtainable, on such a configuration, in the flat- 
spin region. Figures 3 and 4 indicate the following for cases A and B: 


For the slower spin rate (case A) the airplane has oscillatory insta- 
bility for a range of small negative values of C 7 and for all positive 

P 

values of C 7 ; also, for all positive values of C 7 larger than approxi- 

P 

mately 0.05* the motion is divergent. For tne faster spin rate (case B), 
both instability boundaries are moved farther into the positive 

range. In both instances the instability regions are primarily functions 


of C 


P 


This would seem logical inasmuch as 1^ is about 10 times as 


large as 1^ for the airplane considered. (See table I.) Cases A and B 


in figures 3 and 4, respectively, show, then, that an increase in the 
rotational rate could make an airplane having a negative value of 


near zero change from a condition for which the flat spin would be 
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unstable to one in which a stable flat spin would be obtained. This 
effect on the spin stability of increasing the rotational rate has also 
been observed on dynamic models tested in the Langley 20-foot free- 
spinning tunnel. 


Cases C and D when compared with cases A and B show the effect of 
a very small negative value of Ci (-0.01) on the stability boundaries. 

P 

For the slower rotational rate (case C, fig. 5), the region of oscillatory 
instability includes most of the negative Cj range, whereas the diver- 
gence boundary is relatively unaffected by the change in C 7 . Doubling 

P 

the rotational rate (case D, fig. 6) made motions for given values of 
negative Cj stable, whereas they were unstable for case C. 


Cases E to G (table III) were investigated to determine what effect 
a small positive value of C j (0.025) might have on the stability of the 


spin. 


The original values assumed for 


Cmq and c n r 


lead to a stability 


equation which indicated, by inspection, that the system would be unstable. 
(The C2 coefficient in the stability equation (eq. (A21), which contains 

the damping terms was negative.) Accordingly, larger negative values for 
c m q and were as sumed. The stability plots for cases E to G are 


presented as figures 7, 8, and 9, respectively, and indicate that insta- 
bility generally exists for the whole region plotted. As noted by com- 
paring case E with case F , doubling the spin rate r Q and doubling the 

damping in yaw C n ^ caused an extremely small region of stability to 

occur where originally no stability existed. Flattening the slope of the 
pitching-moment curve (decreasing C^ had little effect on the stability 


region obtained for the higher rotational rate. (Compare cases F and G 
In figs. 8 and 9> respectively.) 


On the basis of the results discussed so far, it appears that, for 
an airplane that is loaded predominantly along the fuselage (as in the 
present case), a positive value of Cj in the high-angle-of-attack 

range will tend to prevent stable flat-spinning conditions. For such a 
design, it would be expected that the oscillatory instabilities indicated 
for cases A to G would be evidenced as roll oscillations since the axis 
of the least inertia is the roll axis. 


In order to show the effects of mass distribution on the stability 
of the flat spin, a case (case H) similar to case B was computed except 
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that the moments of inertia about the X- and Y-axes were interchanged. 
Comparison of cases B and H indicates that, although stability existed 
for negative values of when the mass was distributed along the 

fuselage (case B), oscillatory instability is obtained for the whole 
range of positive and negative values of Cj when the mass was dis- 
tributed along the wings. No boundaries are shown in this case (case H) 
because they fell outside the chosen range of Cj . This latter condi- 
tion is undoubtedly attributable to the fact that damping in pitch C m 


was zero. It is probably safe to say that the oscillatory instabilities 
that exist for case H are pitching oscillations, inasmuch as the axis 
of the least inertia is the pitch axis. Spin-model tests have also shown 
that for models loaded predominantly along the wings, the motion is char- 
acterized by pitching oscillations. One conclusion that may be drawn 
from case H is that, for a configuration having its mass predominantly 
along the wings, a value of C nr of not less than approximately -0.10, 

and very little damping in pitch, stable flat-spin conditions will prob- 
ably not be obtainable at the flat-spin attitudes. 


Computed Motions Due to Applicatior of a Forcing Function 

Although the limiting assumptions in deriving equations (A26), (A27) 
and (A28) are such that computation of the motion, after the application 
of any given disturbance, applies only for small changes from the initial 
spin equilibrium position, it was felt that the trends indicated would 
be in the proper sense j therefore, several, computations were made. 

Because of the limitations of this method the time scale has been kept 
down to 20 seconds. 

It appears that the most effective way to influence the spin and to 
bring about recovery (according to refs. 1 and 7) is to obtain a yawing 
moment by applying a moment about an axis which offers the least resist- 
ance to a change in angular velocity (least moment of inertia) . Accord- 
ingly, computed motions were made for applied rolling moments only. 

These computations were arbitrarily made for cases A, B, and E and are 
presented in appendix B. 

The majority of the cases computed apply to case A in table II and 
figure 3 = - 0.10 and r Q = 1-9 radi ins/sec^ . The motions computed 

for the stable case (point 1 in fig. 3) siow that with C nj3 


=0.10 and 
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Cj = - 0.15 when the incremental rolling moment 




originally 


applied to balance was removed, the following occurred: The spin 


rate decreased slightly, 9 became slightly more negative (a decreased 
slightly), and 0 was slightly changed. (See fig. 10.) It should be 

i7~ 


noted that removal of the incremental rolling moment 


in this case 


is equivalent to applying ailerons slightly against the spin, and it 
appeared unusual that ailerons against the spin would slow down the spin 
rate or decrease the angle of attack even slightly for the type of mass 
loading considered. The motion computed for the stable point 2 in fig- 

was taken as -0.10 with C i still equivalent to - 0.15 


lire 5 where 
showed that 


'nr: 


3 

, for the same type of forcing function used for point 1 (that 


is, the equivalent of a small aileron deflection against the spin), the 
spin rate now accelerated somewhat, the angle of attack increased 
slightly (A0 slightly positive), and the angle of sideslip became 
slightly negative. Further, when a rolling -moment coefficient of 0.01, 
which was assumed equivalent to moving ailerons full with the spin, was 
applied for point 2, the combined effects of Cj and ailerons full 

P 

with the spin Indicated a slowing of the spin rate 


a decrease in the angle of attack 


(^- = -0.041 radian/sec 
\dt 

= -0.002 radian/ sec j , and positive angle of sideslip (0 
(See fig. 11). These effects obtained for point 2 where C 


2.9°). 


n, 


P 


was 


equivalent to -0.10 are consistent with the effects that would usually 
be anticipated; whereas the effects obtained for point 1 where CL 

n 0 

was equivalent to 0.10 are inconsistent with the effects that would 
usually be anticipated. 

These differences can be explained on the basis of physical con- 
siderations and also from an examination of equations (Al), (A2), 
and (A3). Physically, it would be anticipated that, for the type of 
loading considered, the cross-couple inertia yawing moment (l x - Iy) pq 

would act in a sense that would slow down the rotational rate when 
ailerons are placed with the spin (refs. 6 and 7) or when a positive 
rolling moment is applied. Similarly, a negative rolling moment 
(ailerons against the spin) acts in an opposite manner. In addition, 
if the spin rate slows down, the angle of attack would be expected to 
decrease because the amount of positive pitching moment supplied by 
the inertia pitching -moment term ^ 1^ - I x jpr would be decreased. 
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Thus, based on Inertia considerations alone, it would be anticipated 
that a positive rolling moment (ailerons with the spin) would slow down 
the spin rate and decrease the angle of attack, with the converse holding 
true for a negative rolling moment. On the other hand, when the aero- 
dynamic yawing-moment characteristics are considered, it would appear 
that, if a spinning airplane had a positive value of C n p, movement of 

ailerons with the spin would cause the inner wing to drop (right wing 
in a right spin) and permit the airplane to acquire a certain amount of 
positive sideslip. This combination of positive sideslip combined with 
positive values of would cause the airplane to yaw into the spin, 

that is, to increase the rotational rate. This effect, then, is just 
opposite to the inertia yawing moment produced by the ailerons. Thus, 
as Cr, increases from zero in the positive direction, there would 

n P 

apparently be some value of positive C n ^ where the aerodynamic and 

inertia effects would be nullified; and for further increases in posi- 
tive C n p, the effects of the ailerons in producing recovery would be 

reversed • 


Examination of the stability equation (eq. (A19) ) and equations (A26), 
(A27), and (A28) for the roll forcing function enables the computation 
of this point of reversal for the cases considered. These equations 
indicate that, for the stable-spin cases considered, the most important 
terms are the coefficients and oi the stability equation 

(eq. (A19)) and the constant term in the roll forcing equation. The 
coefficients and Cg are, of course, both positive for the stable 

cases; and since examination of the roots indicates that the only real 
root is small, A, the equations for Ar and A0 (eqs. (A26) and (A27), 
respectively) for a stable point when a positive is applied can be 

simplified roughly to 
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where K^. is a negative number, K A0 is a positive number, and 0 , 

according to the original assumption, is a small negative angle. Thus, 
the value of M z ^ /or C n j for which application of a rolling moment 

will have no effect on the yawing rate is computed by equating 




( 6 ) 


For case A, this value of C n ^ is 0 . 055 - It is obviously desirable, 

then, that for the application of a positive rolling moment (ailerons with 
the spin) to have the greatest effect in producing recovery for the type 
of loading considered (mass predominantly in the fuselage), C n should 

/ P 

be a large negative number. (it should be noted that, for a loading in 
which the mass is extended primarily along the wings, C n should be 

positive. j It should be pointed out that these values of C n p are 

favorable for instability, which is good in the spin but is not good for 
stability in normal flight. Examination of the roll equations (eqs. (A26), 
(A 27 ), and (A28)) indicates that the wing will always roll in such a 
manner that the sideslip angle acquired is in the same sense as the 
rolling moment ^that is, a positive leads to a positive 0\, 

value of " *- — ' ' ~ 


The 

^np Just computed for which the reversal of aileron effect 


occurs also holds in the oscillatory instability region. It is inter- 
esting to note that incremental negative values of C n p due to the 

fuselage would, in general, be expected when a fuselage nose is pro- 
viding a damping moment in the spin; whereas positive incremental values 
°t *^np would probably be present for cases in which the nose is pro- 
viding a propelling moment. 


Motions for points 3 and 4 in the oscillatory instability and 
divergent regions, respectively, for case A (fig. 3) showed that, for 
the point in the oscillatory instability region (point 3 ), Ar would 
be exceedingly high before any appreciable instability in 0 would be 
obvious; whereas the point chosen in the divergent region (point 4) 
diverged about all three axes almost immediately. (See figs. 12 and 13.) 
In order to note any divergence in 0 for the oscillatory instability 
region, C n p should be close to 0.055> the point where application of 

a rolling moment would have no effect on the yawing rate r or on the 
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angle of attack (90° +0). Although the roots were not computed for this 
case, it is fairly obvious that there would be little change in the roots 
provided remained fixed ^that is, point 3 increased positively to 

^np = 0.055 in fig. 3^. Then, for a point in the oscillatory instability 

region of case A (fig. 3) where CS = 0 and C n R = 0 . 055 ? the motion for 


the case wherein ailerons are placed full against the spin / acting 

\ P 

and = -0.01^ would be approximately as shown in figure 14. As is 

shown, this is a very slow divergence but motions of this nature are often 
seen in the Langley 20-foot free-spinning tunnel. 


A few computations were made for comparison of the motions for the 
various points numbered 1 = -0.15 and C n ^ = 0.1^ for cases A, B, 

and E and show the effect of increasing the rotational rate and of changing 
C 7 from a stabilizing value (-0.10) to a destabilizing one (0.025). The 
P 

most obvious effect for the increased rotation rate is that the response 
to a rolling moment is in the sense anticipated in that a positive rolling 
moment decreases the yawing rate. (See fig. 15 .) This happens because 
the value of C n p for which the effectiveness of the ailerons reverses 

is increased positively from 0.055 for the slower rate of rotation to 
0.155 for the present rate of rotation. The motion for point 1 of case E 
is now unstable (whereas a stable condition existed for case A), and the 
motions plotted in figure l 6 show that the instability in sideslip is 
obvious in a short period of time. 


The effect of applying a positive pitching moment is as follows: 


The change in yawing rate Ar is decreased if is negative unless 

I -r “ T v 

becomes sufficiently positive that 


■'n, 


P 


h : h r n 2 e 

T O O 


m 7 

. h. 


becomes 


negative, p becomes positive provided is negative, and A0 is 


dependent upon various factors and will be positive or negative depending 
upon whether the following expression is positive or negative: 



2 n 2 



% % , 

+ — — — - Q n 
I X X Z ° 


(T) 


Equations (A32) to (A34) show the application of a pitching moment. 
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The effect of applying a negative yawing moment is as follows: The 

change in yawing rate Ar is usually decreased and a negative yawing 
moment has the greatest effect in reducing Ar if and are 

large negative values, (3 generally becomes negative if is nega- 

M X p 

tive , and A9 becomes more negative (a decreases) unless — ti becomes 

Iv 


large enough positively to exceed 


— — r Q 2 . (See eqs. (A29) to (A31) 

Iv 


for application of a yawing moment.) 


Thus, application of a positive pitching moment appears to produce 
an effect similar to the application of a positive rolling moment as 
regards its effect on reducing the spin rate, and a negative yawing moment 
also reduces the spin rate. These effects are consistent with effects 
observed from experimental spin research in the Langley 20-foot free- 
spinning tunnel. 


A few calculations have been made on a high-speed digital computer, 
utilizing six-degree-of-freedom equations of motion, in an attempt to 
check the validity of the results presented in this paper. Although 
some limiting assumptions had to be made in carrying out the six-degree- 
of-freedom studies, the results are considered to indicate the qualitative 
validity of the three-degree-of-freedom results in this paper. 


CONCLUSIONS 


The following conclusions are based on linearized equations of motion 
utilizing only the three moment equations and assuming only flat-spin con- 
ditions : 


1. Contemporary airplane designs (with the moment of inertia about 
the wing axis Iy considerably greater than the moment of inertia about 

the fuselage axis Iy ) having positive values of rolling -moment coeffi- 
cient due to rolling or rolling-moment coefficient due to sideslip will 
probably not have a stable spin near an angle of attack of 90 ° • In addi- 
tion, high rotational rates have a stabilizing effect on flat spins , an 
effect which has been observed during spin-model tests in the Langley 
20-foot free-spinning tunnel. 

2. If the damping in pitch in flat-spin attitudes is zero, it would ' 
not be possible to obtain stable flat spins if the moments of inertia of 
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the airplane are such that the mass is primarily distributed along the 
wings , that is, with I x considerably greater than Iy (opposite to 

the mass distribution of contemporary fighter designs). 


3. The effect of moving ailerons with the spin in producing 
recovery for contemporary fighter designs will be greatest for large 
negative values of yawing-moment coefficient due to sideslip C n ^, an( ^ 

for a certain critical value of positive Cq^ the rolling moment applied 


by moving ailerons with the spin will have no effect on reducing the spin 
rate. These effects also apply to an application of a positive pitching 
moment . 


Langley Research Center, 

National Aeronautics and Space Administration, 
Langley Field, Va., February 25, 1959* 
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APPENDIX A 

DERIVATIONS OF EQUATIONS 


Equations for Determining Stability Boundary 

The equations of motion corresponding to the rolling, pitching, and 
yawing degrees of freedom are, respectively. 


Roll 


I xP = 

( Z Y " I z) ( l r + M X p P + M XpP + ^X 

(Al) 

I Y q = 

( z z " z x)p r + m y 0 0 + % q <i + 

(A2) 

: z f " ( 

;i x - I Y )pq + M z ^r + M Zp P + AM Z 

(A3) 


Expressing the angular velocities about the body axes in terms of 
the (Euler's) attitude angles gives 

• • 

P = 0 - sin 9 (A4) 

q = 0 cos 0 + V cos 0 sin 0 (A5) 

r = -0 sin 0 + ijr cos 0 cos 0 (a6) 

(See appendix A in ref. 8, for example.) The following assumptions 
were applied: 


0 = 0 O + A0 
V = r Q + Ar 

P = 0 = P D + Ap (P 0 = °) 

sin Z = Z- 
cos Z = 1 

Rewriting equations (A4), (A5), and (A 6) with the aforementioned 
assumptions included and with second-order terms neglected gives , 
respectively, 
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p = A|3 - r Q e o - r Q A0 - e Q Ar 

(A7) 

q = A0 + r Q AP 

(A8) 

r = r Q + Ar 

(A9) 

which can be expressed as 


11 

fc: 

1 

O 

fc. 

1 

CD 

O 

(A10) 

q = A9 + r Q AP 

(All) 

r = Ar 

(A12) 

or as 


pr = r D Ap - 2r o 0 o Ar - r Q 2 Ae - r Q 2 0 o 

(A13) 

pq = -r Q e o A9 - r o 2 0 o AP 

(Aik) 

qr = r Q A0 + r ( 2 AP 

(A15) 


Substituting equations (A7) to (A15) into equations (Al) to (A3) 
and neglecting second-order terms gives the three moment equations in 
the following form: Roll is given as 



(A16) 



where 


and 


1 ? 


M- 


lX P 


Iv r ° 9 ° I 




cancel each other for balance at 3=0° 

L X x X 

for the initial steady spin. Pitch is given as 


Pitch = Ap 




A9 


M y 

D 2 - D + 


I Z “ Z X 2 
7 *o 


M, 


Ar 


2r 0 8 0 


I Z * I X 


My p I 7 - Iy M Yq 

-r o 2 0o + ^ 9 o + 


” 1.0 


(A17) 


My 

where — is considered to be equivalent to zero for the initial steady 
X Y 

/ I - I ra Y 

spin and ( -r o 2 0 o -? !— — £ 




cancel each other for balance 


for the steady-spin state. Yaw is given as 


Yaw = Ap 


X x - % r 2 C % 

+ A6 

■ - 1 *r 9 T) 

l 4- Ar 

n _ 1 

I Z ° ° I Z J 

t r 0 W 0 U 

_ z 

T /AI 

I •% 

z z_ 


— - r + — - 

I z ° I z 


(A18) 


% _ ^2 


where r n and 

h x z 

steady spin. 


cancel each other for balance for the initial 
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By placing equations (Al6), (A17), and (Al8) in determinate form 
and solving for the characteristic equation, an equation of the following 
form is obtained: 

CD 5 + + C + C, D 2 + C D + C, = 0 (A19) 

12 3^5d 

where the coefficients of the characteristic equation are given as 



4 V M Y M v My H, M y 

i + 2a . 2i + -h 


h h 


T Z I X 


I Z I Y I Z 


(A20) 

(A21) 
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Equations for Solution of Motion After Application of 
Various Forcing Functions 

In order to solve for the motion afte:' various forcing functions 
were applied, the five roots of the characteristic equations were deter- 
mined as wore the differential equations for the various forcing func- 
tions. The solution for any of the motions A9, Ar, or Ap is then 
determined from 


A9, Ar, or AP 


-(d) 

: ? (d) 


where f(d) is the numerator of the equation for the motion after various 
forcing functions are applied and F(d) is the characteristic equation. 
The forcing functions employed were step Inputs for which Heaviside's 
expansion could be employed to obtain a solution 


A0, 


Ar , or 


Ap = 


f(o) 

F(0) 


+ 



fQO p At 
XF'(^) 


where the various values of A are the five roots of the characteristic 
equation and F' indicates the first derivative of the characteristic 
equation. In order to solve this equation for the complex roots, the 
method explained in reference 9 was employed. The equations for the 
motion after application of various forcing functions and the removal 

of fhe term originally inserted intc equation (Al6) for balance at 

I X 

zero sideslip are as follows: 


M- 


Ar = 


Xp 


r o®o 





w 

+ I v I„ D 


+ 
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AQ 



r 0 
o o 



Iz ' Ix V g + 

h I 




2r o 5e o 2 




+ 2r Q 0 o 




(A27) 



2 2 

2r o 9 o 




Mz rAz ~ 1 

I Z \ I Y 


X 




(A28) 


Application of a rolling moment .- After the application of a rolling 
moment the equations for Ar > A0 9 and AP are the same as those given 
in equations (A2 6 ), (A27), and (A 28 ), respectively } except that the 

M Xp AM y 

— — r Q 0 o term is replaced by - — 

Z X X X 

Application of a yawing moment .- After the application of a yawing 
moment the equations for Ar, A0, and Zip are given, respectively, as 



2k 


M*y 

Ar = — 


My M v \ /M v M v M y M v 
d 4 - f A + JiV + I i Y( i x p Y e 


L 


X X X Y 


I X X Y I X I Y 


r2 f ± i z - i x i Y - ^ 2 + /% % + My e % 

° \ I Y X X // \ X X X Y % % 


- r, 


My 

2 1( l 


. 2 %l +r2 ^% + “ 

O T / O T T T J. 

1 X / x X X Y 


0 X Y “ X Z 2 I Y " I Z I Z " I X 4 
i _ j j r Q 

i X X Y 


Y P I v - I. 
Y 


% ° 


M Y T T O M Y 

x p I z ~ I x „ 2 ^ Y e X P 


r„ + 

O T T 

X Y x X 


1 

FldJ 


(A29) 


My 

A© = / 0 o r o 

X Z 


/ T T \ /My M v M-v 

I Y / Vy I X X X 


*p , \ x z - 4 X] 


D + 


I Z “ I x/ I Y T Z 2 Mx 3\ Mx p| 1 


r + 

% ° x x 


X Y Z X 


F(d) 


(A30) 


^Z 

AP = ~ 0 O 

i Z 


3 I V, "M„2 


D- + [- IZ I IX r 2 - 2 Iz ' Ix I r 2 + 

I v ' I v / I Iy ° Iy I X ° 


D “ | f— 1 + 

\ X Y - L X 


% %■ *1, 


I x I Y Iy 


My T T M. My 

10 ) D + __P h ~ x X r 2 + *p Y 0 


Iy 


F(d) 


(A31) 
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Application of a pitching moment .- After the application of a 
pitching moment y the equations for Ar, A9, and Aj3 are given , respec- 
tively, as 
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APPENDIX B 


EQUATIONS FOR MOTION OBTAINED WFEN VARIOUS FORCING 
FUNCTIONS ARE APPLIED IN ROLL 


Case A 

In the determination of points 1 to 4 in figure 3 for case A, the 

AM Y 

_ term is removed from equation (Al6) wl ich then allows Mjr ac ^' 

I x ** 

Thus, Ar, 00, and Ap can be expressed, respectively, as follows: 
Point 1 in figure 3 - - 

Ar = -0.03903 + 0.039064e -0 '° 2 ^ t + 0.000196e _ 0 , 017 l 84 cos(l .696t + 
6.1665) + 0.00005l4e‘°‘ 150l|5t cos(3.32at + 3-6249) 


A0 


= -0.00252 + 0 . 002521e -0 • °267 t + o.0004e -0-017l8t cos(l.696t + 3-09) + 


0 . 0000156 e _ 0 , 15 ° 45 t cos( 3 . 328 t + 2-904 1 

AP = -0.005466 + 0.010T36e'°- 0267t + 0 .000628e - 0 ’ 017 l 8 t cos(l.696t + 

4.606) + 0.001126e“° ,15045t cos(3-328t + 4.647) 

Point 2 in figure 3 - - 

Ar = 0.128656 - 0.12879'Se'°‘ ° 2925t 


A8 = 0.00946 - 0. 0094 7e ‘°* 029254 
AP = -0.004996 - 0.0002 71e - ° ' 02 925 4 

It should be noted that the oscillatory terms were not computed for this 
case in that they were assumed to be of negligible importance for the 
time range considered-. 
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If C l = 0.01 is applied, the result is 

Ar = -I.3291 + 1 • 3305e _0 ’ 02 925 1 

A0 = -0.09772 + 0.0978e“°- O2925t 

Ap = 0.051648 + 0. OO276 9 e ' 92929 ^ 

Point 3 in figure 3 .- 

Ar = 0.5764 - 0.5878e-°-° 95t + 0.0058e-°- 17t cos(l.32t + 1 
0.00092e 9 ' 9129t cos(2.05t + O.729) 

AS - 0.0424 - 0.0432e-°-° 35t + 0.00026e~°- 17t cos(l.32t + ( 
0.0001e°' 0125t cos(2.05t + 1.84) 

AP = -0.02346 - 0.005e-°-°55t + o.029e" 0 - 17t cos(l.32t + 6. 
0.0004e 0,0129t cos(2.05t + 5.28) 

If C l = -0.01 is applied, Ap can be written as 

AP = -0.2438 - 0.0558e-°- 059 t + 0 .1122e-°- 17t cos(l.32t + 6 

0.0027e°' 0125t cos(2.05t + 5.14) 

Point 4 in figure 3 .- 

Ar = -0.05515 + 0.00028e 2 ' 798t - 0.00022e _3 ' 13t + 0.05013e 
0.00002e°- 005t cos(2.043t + 4.94) 


•37) + 


•99) + 


18) + 


• 19) + 


0.03t 
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A9 = -0.00405 + 0.00007e 2 ‘^ t - 0.0000 6« _5 ' 15t + 0.00368e °'° 5t + 

0.00013e 0 ' 00 ^ t cos(2.043t + 1.621) 

AP = 0.00593 - 0.00309e 2- ^ t - 0.00275e“^’ 15t + 0.00044e °'° 5t + 

0.00003e°' 005t cos(2.043t + 5-554) 

If Cj = -0.01 is applied, AP can be written as 

AP = -0.05651 + 0.002l66e"°' 02 ^ t +0.0348e~°* 01 ^ l8t cos(l.696t + 6.19) + 
0.03876e" 0,15 ° 45t cos(3.328t + 6.26) 


Case B - Point 1 in i 'igure 4 


In the determination of point 1 in figure b for case B, the ^ 

X 

term is removed from equation (Al6) to allow to act. Thus, 

Ar = 0.0597 - 0.0596e-°- 02 ^ t + 0 .0000274e"°‘° 565t cos (3 - 133t + 4.41) + 


0.000228e -0,22 ^ t cos(3.65t + 1.60) 


A0 = 0.003 - 0.0029e"°‘° 2 ^ t + 0.00092e 


-C.0363t cos ^ <155t + 1>l8 ) + 


0.00078e _0 ‘ 225t cos(3.65t + O.78) 

= -0.0066 + 0.000056e“°' 02 ^^ t + O.OOOS€ -0 *°^^ t cos(3.133t + 2.63) + 


0.0055e"°' 225t cos(3-65t + 6.07) 

If Cj = 0.01 is applied, Ar, A0 , and Ap (with oscillatory 

terms being of small importance in determining A0 and final trim P) 
are given, respectively, as 
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Ar = -0.3T82 + 0.3776e _0 * 02T5t + O.OOl648e _ °'° 565t cos(3.133t + 1.27) + 
O.OOl442e“°' 225t cos(3.65t + 4.73) 

A0 = -0.016 + 0.0l6e" 0,02 ^ 5t 
Ap = 0.035 + 0.0003e"°- 02 ^ t 

Case E - Point 1 in Figure 7 

AMx 

In the determination of point 1 in figure 7 for case E, the 

I x 

term is removed from equation (Al6) to allow to act. Thus, 

Ar = 0.002707 - 0.002703e _0 -° 92t + 0.0000203e"°’ 005t cos(2.0l8t + 4.71) + 
0.000092e 0 *° 59t cos(3.15t + 1.62) 

A9 = 0.000204 - 0.000208e“°-° 92t + 0.0001e“ 0 '°° 5t cos(2.0l8t + 1.53) + 
0.000067e°-° 59t cos(3.15t + 4.71) 

AP = 0.00131 - 0.000055e“ 0,092t + 0.000049e -0 ' 00 ^ t cos(2.0l8t + 0.04) + 

0.00131e 0,059t cos(3.15t + 3-17) 

If Cj = -0.01 is applied, the result is 

Ar = -0.111864 + 0.111732e“ 0,092t + 0.00084e“°’ 005t cos(2.0l8t + 1-57) + 

0.0038e 0,059t 


cos(3.15t + 4.76) 
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A0 = -0.00839 + 0.00856e" 0,092t + 0.00472e~ 0 *°° 5t cos(2.0l8t + 4.67) + 
0.00276e 0,059t cos(3.15t + 1.67) 

A3 = -0.054337 + 0.002268e*°'° 92t + 0.00202e' 0,0 ° 5t cos(2.0l8t + 3-18) + 
0.054256e°*° 59t cos(3.15t + 0.01) 
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TABLE I.- MASS AND DIMENSIONAL CHARACTERISTICS OF AIRPLANE 

l x , slug-ft 2 10,000 

l y, slug-ft 110,000 

l z , slug-ft 2 115,000 

S, sq ft 385-33 

b, ft 35.7 

Test altitude, ft 30,000 



TABLE II.- CONDITIONS INVESTIGATED 




























































TABLE III.- COEFFICIENT'S OF CHARACTERISTIC EQUATION FOR 
THE CASES INVESTIGATED (SEE TABLE II) 


Coefficient 

Cane A 

Cane B 

Cane C 

c i 

1 

1 

1 

C 2 

0.35 

0.35 

0.06l 

C 3 

0.371C np - 53C Zp + 6.022 

0.371C n(3 - 33C lp + 13-506 

0.3YlC np - 33C 2p + 6.029 

c 4 

0.125C n „ - 1.484C 7 + 1.620 

P t p 

0.125C np - 1.484C l} + 3-813 

0 .012C n _ - 1.484C 7 + 0.508 

n P £ p 

C 5 

0.260C rt{B - 219. W 2 + 7-314 

0 . 260C - 540.826c 1 + 50.356 

0.260C n - 219. 95C z p + 7-301 

C6 

-0.339C np - 6.157Cjp +0.229 

-1.130C- - 14.673c, + 1.468 
n P " l p 

; -0.032C n - 5.247C. + 0.212 



Coefficient 

Cane D 

Cane E 

Case F 


c i 

1 

1 

1 


c o 

0.061 

0.02 

0.113 


c 3 

• 0.371C np - 33C 1 + 22.090 

0.37ic n -53C, - 6.006 

P ; P 

0.371C nR - 33C 1 + 22.073 
P P 


C 4 

0.012C n _ - 1.484C 7 + 1.112 

n p £ p 

-0.028C n p - 5 • 353C 1 ( + 0.239 

-0 . 028C n - 10. 282C 7 +2.980 

n P £ P 


c 5 

o.26oc np - 702.25 C ip + 12-554 

0.265C np - 220.044C , + 7.447 

0.265C np - 774.94c I + 103.313 


c 6 

-0.139C rip - 20.67C lp + 0.357 

0.085C n{3 - 20.23C Z} + 0.685 

0.408C n(J - 142. 305C 1 + 19.316 


Coefficient 

Case G 

Case H 

C 1 

1 

1 

C 2 

0.113 

0.057 

C 3 

0.371Cn p - 53C lp + 21.738 

0- '>8^C n - ^-Siec, + 21.957 

p p 

C4 

-0.028C np - 10 . 282C 1 + 2.944 

0. )llC np - 0.13 1 >C i(j + 1.113 

c 9 

0.132Cnp - 757- 7l6C i p + 102.913 

4-5 ’9Cn p - 59-2l8C Ip + II6.368 

c 6 

Cj4l9C nn - 1 38. 99C i „ + 18.861 

0. >28c n - 1.656c, + 3.355 


P P 

P P 
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(b) 0 and 0 = 0. 



(c) 0 and = 0. 

Figure 1.- Body system of axes and related angles. 
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+M Y,o 



Figure 2.- Sketch indicating relative nature of various terms in 

pitching -moment equation. 


Figure 3.- Stability diagram for case A 





























Figure 10.- Motion for point 1 of case A (fig. 3 ) with C 7 acting (no rolling -moment coef 





Figure 14.- Motion for point 3 of case A (fig. 3 ) after application of a rolling-moment coef 

ficient of -0.01. (Note the time scales.) 
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